The nature of the magnetic order in Ca3Co206 
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We present a detailed powder and single crystal neutron diffraction study of the spin chain com- 
pound Ca3Co2 06. Below 25 K, the system orders magnetically with a modulated partially disordered 
antiferromagnetic structure. We give a description of the magnetic interactions in the system which 
is consistent with this magnetic structure. Our study also reveals that the long-range magnetic 
order co-exists with a shorter range order with a correlation length scale of '^ 180 A in the ah plane. 
Remarkably, on cooling, the volume of material exhibiting short range order increases at the expense 
of the long-range order. 

PACS numbers: 75.30.Gw, 75.25.+Z, 75.30.Fv, 75.50.Ee 
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The appearance of plateaux in the magnetization 
curves of low dimensional quantum spin systems (e.g. 
NH4CUCI3 [H and SrCu2(B03)2 [2,]) has generated con- 
siderable attention from both an experimental and a the- 
oretical point of view. In this respect the spin chain sys- 
tem Ca3Co206 is very interesting because at low tem- 
peratures several steps, equally spaced in magnetic field, 
appear in the magnetization [j, y, l5[ , a behavior reminis- 
cent of quantum tunneling of magnetization in molecular 
magnets [6|. The origin of this intriguing phenomenon 
is still an open question [71, Isf and the magnetic proper- 
ties of Ca3Co206 have been intensely studied in t he p ast 
decade using many techniques, including x -ray [9|, [lOl, UjJ 
and neutron 12|, [isl, [ij] scattering, NMR [l5| , calorime- 
try [l6| , and magnetometry [3|, [j, |5| . 

Ca3Co206 is a rare example of a material where ferro- 
magnetic (FM) ID Ising spin chains are coupled through 
a much weaker antiferromagnetic (AF) exchange on a 
triangular lattice. The spin chains in Ca3Co206 are 
made up of alternating face-sharing octahedral (Col) 
and trigonal prismatic (Coll) CoOe polyhedra, running 
along the c axis and arranged in a triangular lattice in 
the ab plane |l7|. The different Co environments leave 
the Co^+ ions on the Col sites in a low-spin {S=0) 
state, and those on the Coll sites in the high-spin {S=2) 
state [9I, UjJ, ll5| . Crystalline electric fields also lead to a 
very strong anisotropy with the moments preferentially 
aligned along the c axis [3|, |^ . The zero-field magnetic 
structure of Ca3Co206 has yet to be described unam- 
biguously. Most theoretical descriptions of the mag- 
netic structure center around the two stable configura- 
tions for ID Ising chains coupled antiferromagnetically 
on a 2D triangular lattice, the ferrimagnetic structure 
(M, M, — M) and the partially disordered antiferromag- 
netic (PDA) structure (M, -M,0) [18] (the labels in 
parenthesis indicate the relative magnetizations on the 



triangular lattice). The step-like magnetization has been 
described using the Ising model on a triangular lattice [3] . 
There are, however, experimental data which suggest a 
more complex magnetic structure for Ca3Co206. A pro- 
nounced drop in the intensity of the magnetic peaks on 
cooling below 18 K in zero field has been observed in pow- 
der [12I, ll3| and single crystal diffraction studi es Il4^ . Re- 
cent resonant x-ray scattering (RXS) studies [lO[ showed 
a small (~ 0.01 A~^) incommensuration in the magnetic 
reflections associated with a long-period modulation of 
the magnetic structure along the c axis. 

In this work we address three fundamental questions 
that are key to understanding the physics, including 
the step-like behavior of the magnetization at low T, of 
Ca3Co206. 1) What is the true magnetic ground state of 
this material? 2) What is the physical origin of the ob- 
served modulation in the magnetic structure? 3) What 
produces the low T reduction in the magnetic intensity 
seen in previous RXS and neutron scattering studies? 
To this end, we have carried out a careful investigation 
of Ca3Co206 using neutron diffraction measurements as 
a function of T. The results reported here reveal that the 
actual magnetic order is neither a simple PDA nor ferri- 
magnetic, but corresponds to a longitudinal sinusoidally 
modulated structure with a very long periodicity. We 
present a description of the magnetic interactions in this 
system, that can explain the observed modulation and 
emphasizes the need to consider the 3D character of the 
magnetic couphng in Ca3Co206. Finally, our data reveal 
that the reduction in the magnetic intensity of the neu- 
tron diffraction peaks is due to a coexistence of long and 
short range magnetic order. The results presented below 
provide a new insight into the nature of magnetic order 
in this geometrically frustrated material and call for a 
comprehensive revision of the theoretical models used to 
describe the magnetic behavior in Ca3Co206. 
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FIG. 1: Single crystal neutron diffraction data for Ca3Co2 06 
taken on the time-of-fiight SXD instrument at ISIS ^19]. Neu- 
tron intensity patterns were recorded in the (/i/cO) scattering 
plane at three different temperatures (a) T=30 K, (b) 20 K 
and 2.2 K (not shown), above and below the magnetic or- 
dering temperature, Tat =25 K. A refinement of the data [20] 
(;^2^g j^g^ i^p.2=18.84, i?^^2=24.53) revealed that no struc- 
tural transitions occur on cooling from 30 K down to 2 K; 
there are no significant variations in the lattice reflections or 
the atomic parameters, which are in good agreement with 
those in the literature [12, [l3 • 



Single crystals of Ca3Co206 were grown in a KCO3 flux 
using Ca3Co206 powder synthesized via a solid state re- 
action y. A single crystal 15 x 2 x Imm^ with the longest 
direction parallel to the c axis was used for the neutron 
diffraction experiments. The high quality of the crys- 
tals was confirmed by x-ray diffraction, energy dispersive 
x-ray, magnetization, and specific heat measurements. 

Neutron single crystal diffraction experiments were 
performed on the SXD time-of-flight instrument at the 
ISIS-RAL, UK. SXD uses the white beam Laue tech- 
nique, a stationary crystal, and a large area position- 
sensitive detectors covering a solid angle of ~ 27r sr [22| , 
allowing a quick data collection over a large area of re- 
ciprocal space. Neutron powder diffraction experiments 
were conducted on the GEM instrument (ISIS) and were 
used to follow the T evolution of the magnetic structure. 
The 4-circle TriGS thermal neutron single crystal diffrac- 
tometer at the SINQ facility of the PSI, Switzerland was 
used for precise measurements as a function of T of the 
position, intensity, and width, of a selection of magnetic 
and nuclear peaks. With a graphite (002) monochroma- 
tor providing a wavelength of 2.32 A, the TriGS resolution 
Ad/d was about 0.5%. 
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FIG. 2: (color online). A selection of the results of the pow- 
der measurements on the time-of-flight GEM diffractometer 
at ISIS, (a) Neutron intensity of the (110) nuclear Bragg peak 
in Ga3Go2 06 at different temperatures; no variations in in- 
tensity are observed except for an increase of the background 
at T=30 K due to a paramagnetic contribution, (b) Neutron 
diffraction pattern at 18 K (• j and refinement shown for two 
values of the propagation vector (red dotted Une, kz = 1.00) 
and (blue dashed line, kz — 1.01) [231 ]. 



The results of the SXD measurements are shown in 
Fig. [H Magnetic reflections appearing below Tat =25 
K can be indexed by the propagation vector k= (0,0,1) 
with respect to the he xag onal setting of the rhombohe- 



dral space group R3c [23|]. The only apparent difference 
between the data collected at 20 and 2.2 K is an ~20% 
drop in the intensity of the magnetic peaks at low T, as 
previously reported 12|, ll3|, |lj] . All the (hkl) magnetic 



peaks have even /, which supports the idea that the intra- 
chain ordering is ferromagnetic. The systematic absence 
of the (00/) class of reflections confirms the Ising nature 
of the magnetic system with the Go magnetic moments 
aligned parallel to the c axis. Previous neutron diffrac- 
tion studies have suggested Ga3Go206 adopts a ferrimag- 
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netic structure 



|. This assumption can be rejected 
as this would lead to additional magnetic contributions 
at k=0, and especially in the (110) and (300) peaks, that 
are not observed experimentally in either single crystal or 
powder neutron measurements (see Fig. |2^). Our SXD 
single crystal data is consistent with a simple PDA struc- 
ture. The refined value at T=20 K of the magnetic mo- 
ment on the high-spin {S=2) Coll ion is 5.0 ±0.1 fiB [24 1 
and includes a sizeable orbital moment (L~l). Magnetic 
dichroism measurements on Ga3Go206 gave similar val- 
ues for the spin and orbital moments on the GoII ion [9| . 



GEM powder diffraction data, of slightly higher res- 
olution, reveal that the magnetic propagation vector 
is in fact incommensurate with the crystal structure, 
k=(0,0,/c2) kz ~1.01, in excellent agreement with recent 
results obtained with magnetic x-ray scattering [lO|. 
This is illustrated in Fig. [2)3, where a commensurate 
value for kz fails to reproduce the correct position of 
the magnetic Bragg peaks. Symmetry analysis using 
representation theory shows that the magnetic represen- 
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FIG. 3: (color online), (a) Schematic of the proposed incom- 
mensurate magnetic structure for 30 unit cells along the c axis 
with kz increased to 1.05 to emphasize the modulation of the 
moments. The magnitude and direction at sites in three ad- 
jacent chains are shown, (b) Magnetic exchange interactions 
in Ca3Co2 06 Ji (black lines) FM; J2 (red dashed lines) AF; 
J3 (blue dotted lines) AF. Coll sites at (0,0, |) and (0,0, |) 
are labeled 1 and 2. 
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FIG. 4: (color online). Transverse scans in the ab plane 
through the (100) magnetic Bragg peak at (a) T=5 K and 
(b) r=21 K. The data were collected on the 4-circle TriCS 
diffractometer at PSI. Each peak is fitted (full line) by the sum 
of a broad and a narrow Gaussian function, together with a 
flat background, (c) FWHM and (d) integrated intensity of 
the broad (A) and narrow (■) Gaussian components of the 
(100) peak as a function of T. The total integrated intensity 
(•) is also shown. 



tation F for a magnetic ion on site 6a (0,0, |) is de- 
composed into three irreducible representations: F = 
Fi + F2 + 2F3 [25|. Only the symmetry-adapted mode 
belonging to Fi allows a fit to the data. The magnetic 
phase transition, involving a single irreducible represen- 
tation, is therefore consistent with the Landau theory 
of second-order transitions. For this mode the moments 
lie along the c axis and their amplitude can be calcu- 
lated using Mi(Rl) = Mcos (27rk • Rl) and M2(Rl) = 
—Mcos (27r (k • Rl + "^^2^)) where the subscript 1 and 2 
refer to symmetry related sites 1:(0,0, |) and 2:(0,0, |), 
and Rl = {Rx-,Ry,Rz) is the translation vector with 
respect to the zeroth-cell including R-centering transla- 
tions. The resulting magnetic arrangement with an ex- 
tended longitudinal modulation, (~1000 A or 200 Coll 
ions) is illustrated in Fig. [3^. If the periodicity is truly 
incommensurate with the nuclear lattice, (kz is not a 
fractional value), any value of moment between -M and 
+M is found somewhere in the lattice. At specific lattice 
points, the magnetic configuration within nearest neigh- 
bor sites in adjacent chains of the triangular lattice is 
exactly the simple PDA structure with (M, — M, 0) con- 
figuration. At other lattice points, the configuration is 
exactly (M, -M/2, -M/2). Other intermediate situa- 
tions are also found. This type of structure is usually 
induced by competing exchange interactions in the pres- 
ence of strong axial anisotropy, stabilizing an arrange- 
ment in which the ordered component of the magnetic 
moment fluctuates along the propagation direction. 

Ca3Co206 is often described in terms of Ising chains 
on a triangular lattice. Within this framework, the AF 
in-plane correlations and the FM correlations along the 



c axis are completely decoupled. Clearly, the magnetic 
energy of such a lattice is not lowered by introducing a 
small incommensurability, so the approximation to a 2D 
lattice must be revised. The Co ions in adjacent chains in 
Ca3Co206 are shifted out of plane by | or | of c |17|]. A 
sketch showing the interactions within a single triangular 
unit in the structure and omitting the oxygen atoms for 
simplicity, is shown in Fig. [Sb. Within each chain, direct 
Co- Co overlap dominates and leads to a strongly ferro- 
magnetic interaction Ji [26|. The inter-chain coupling 
is more complex. There are two inequivalent AF super- 
superexchange interactions involving Co- 0-0- Co paths 
labeled as J2 and J3, that follow helical paths and con- 
nect the Co sites in adjacent chains. For J2, the overlap 
of the O 2p orbitals is very small making J3 the dom- 
inant term after Ji [26|. This helical exchange path- 
way introduces an AF inter-chain coupling between the 
nearest neighbor Coll ions that competes with the intra- 
chain FM interaction producing the observed modulation 
along the c axis. This type of competition leads nat- 
urally to non-colinear helicoidal arrangement for weak 
single-ion anisotropy, or a longitudinal sinusoidally mod- 
ulated structure in the presence of large axial anisotropy 
as described here. We note that the present structure 
contains extended regions with very weak static mag- 
netic moments, a situation that favors the development 
of defects or dislocations along the chains. 

Finally, we address the origin of the reduction in in- 
tensity of the magnetic peaks. Figures [4^ and [4)3 show 
transverse scans in the ah plane through the (100) mag- 



netic Bragg peak as measured on the TriCS diffractome- 
ter at T=5 K and T=21 K. The peaks are fitted us- 
ing two Gaussian functions with significantly different 
widths, the narrower of which corresponds closely to the 
instrumental resolution. The T dependence of the peak 
widths are shown in Fig. [It . Above 21 K, the peak can- 
not be fitted reliably with the two components, therefore 
the data are truncated at this point. At T=21 K, the 
ratio of the intensities of broad and narrow components 
is 0.2; this ratio increases to 0.7 at low T. A similar 
behavior has been observed for several other magnetic 
peaks, including the (200), (112) and (102) and in differ- 
ence plots of the GEM powder diffraction data collected 
at different temperatures. The loss of integrated inten- 
sity in the narrow component is mirrored by an increase 
in intensity of the broader feature as well as in the overall 
background. In other words, the anomalous reduction in 
the intensity of the magnetic peaks is due to the onset 
of a short-range magnetic ordering. From the reflection 
widths measured at 5 K, the estimated correlation length 
for the short range order is 180 A in the ah plane [27|. 
This is a clear sign of the increasing instability of the 
longer-range magnetic order as the temperature is re- 
duced. One can speculate that the onset of shorter-range 
correlations originates from defects in the propagation of 
the long-range magnetic structure presented earlier, with 
stacking faults developing along the c direction. These 
results are clearly different from the case of CasCoRhOe 
where short range magnetic order is observed above and 
below the long range ordering temperature of 100 K, but 
there is little or no T dependence to the diffuse compo- 
nent below 100 K 0. 

In summary, neutron diffraction measurements show 
that the magnetic structure in Ca3Co206 is a complex 
modulated PDA structure running along the c axis. We 
propose a description of the magnetic interactions in 
terms of helical path couplings. This 3D description can 
explain the observed modulation along the chains. Our 
data reveal the appearance, at low T, of a short range 
magnetic structure that coexists with the long range or- 
der. We suggest that this phase coexistence is the origin 
of the drop in the intensity of the antiferromagnetic re- 
flections at low T observed in previous experiments. 
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